NACA TN 3555 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3555 


A METHOD FOR CALCULATING THE CONTOUR OF BODIES OF 

REVOLUTION WITH A PRESCRIBED PRESSURE GRADIENT AT 

SUPERSONIC SPEED WITH EXPERIMENTAL VERIFICATION 

By Paige B. Burbank 

Langley Aeronautical Laboratory 
Langley Field, Va. 



W ashington 
March 1956 





NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 5555 


A METHOD FOR CALCULATING THE CONTOUR OF BODIES OF 
REVOLUTION WITH A PRESCRIBED PRESSURE GRADIENT AT 
SUPERSONIC SPEED WITH EXPERIMENTAL VERIFICATION 
By Paige B. Burbank 

SUMMARY 


A method based upon the theory of characteristics has been developed 
to compute the contour of a body for a prescribed pressure distribution. 

An indication of the accuracy of the method was obtained by using the 
pressure distribution which had been determined by the conventional method 
of characteristics to compute the ordinates near the nose of a body. The 
agreement between the original and computed ordinates was good. A model 
was constructed with an initial included cone angle of 50^ and a contour 
to give a linear pressure distribution for a free-stream Mach number 
of 5 . 13 . Good agreement between the experimental and computed pressure 
distributions was obtained when the model was tested at free-stream Mach 
numbers of 3*05 and 3*13. 


INTRODUCTION 


In the past, no method that permits computation of a body of revo- 
lution with a particular pressiire distribution has been available. The 
design of such a body would be useful in the investigation of many aero- 
dynamic problems such as those associated with body-wing interaction, 
inlets, and boundary -layer transition, separation, and shock interaction. 
In this paper, a method based upon the theory of characteristics has 
been developed to compute the contour of a body for a prescribed pres- 
sure distribution. This method has been used to compute the contour of 
a body which has a constant pressure gradient. This pressure distribu- 
tion was chosen since the theoretical approach to the problem of tran- 
sition indicates that it would be desirable to have, for certain experi- 
mental studies, a constant pressure gradient. The available experimental 
data on bodies of revolution is restricted to cones which have a zero 
pressure gradient and to bodies which have a varying pressure gradient. 

A scale model has been constructed and the experimental pressure distri- 
bution checked with theory. 
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SYMB0I5 


Cp pressure coefficient 

I image 

^ _ sin (j. tan p sin 9 
cos(0 + p) 

L length of body 

distance from nose to point of r ^^^ 


m = 


M 

0 


P 


Po 

5>o' 


0. = 


sin [g tan ^ sin Q 
cos(0 - 1^) 

Mach number 

object 

static pressure 
stagnation pressure 

stagnation pressure behind bow shock 



R 

V 

Vz 

w 


radius 

Reynolds mmiber 
velocity 

limiting velocity 

limiting velocity ratio, V/V^ 


x,y 


coordinates of flow field 


a = 


9c ^ 6a 


3 = e - 6 
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6 = 0g + dGg^h 


e 

shock angle 

7 

ratio of specific heats^ 1.405 


Ai = ^ = tan(e + n) 
Aji = ^ = tan(e - n) 


n 

conical flow angle 

0 

stream direction 

0 

meridian angle about axis of symmetry 

1^ 

Mach angle 


Subscripts: 


b 

body 

c 

cone 

m 

mean value; when preceded by a 2^ a second mean value is 
indicated 

max 

maximum 

n 

any point 

N 

normal to shock 

t 

theory 

00 

free-stream conditions 

a y b y C y • • • 

points in the characteristic net 


METHOD OF DESIGN 


A sharp-nose body of revolution vlth a prescribed pressure gradient 
may be designed in the following way: The nose can be considered conical 

for a short distance^ and the values of the flow properties - |i the Mach 
angle and 0 the stream direction - along the characteristic line af 
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(see fig. l) can be determined for any conical flow angle ti from refer- 
ence 1. If Cp is assumed to be linear with x (see fig. 2), the slope 
of the body at any point is determined by the method of characteristics 
as outlined in reference 2. For a characteristic of Ajj (the entropy 

terms being disregarded), the basic formulas are as follows; 


^ = tan(e - M-) 


( 1 ) 


^ + de tan p - m ^ = 0 
V y 


( 2 ) 


If the limiting velocity ratio 

Tr 2 


'I _ 


V2 


= 1 + 




7 - 1 


(5) 


is introduced and equations (l) and (2) are applied to the determination 
of the flow properties at a point c on a specific body (figs. 1 and 2) 
with the following relations: 

dW = Wc - Wb 


<ieb,c = 0c - 6b 

^^b,c = ^c - ^b 

where 0 is measured in radians and equation (2) becomes 




”^b *^^b,c ^c - ^b 


•^0b,c ” 


yb 


Wv 


tan p-|^ 


(5) 


where 


m = 


sin 6 sin p tan p 
cos(0 - p) 


(6) 


the 0(3 corresponding to each value of d0-b^c is 


0c - ^®b,c ®h 


(7) 
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The value of and Wq to be used in equation ( 5 ) is determined 

as follows: Draw a straight line from point a to an intersection of 

the 7\jj characteristic originating at point b. For an approximation, 
the angle of the line to the horizontal is slightly less than the slope 
of the conical surface leading to point a. This intersection is brack- 
eted by several values of dx-|^ ^ and their corresponding values of 

are determined from the Cp variation given in figure 2 and the following 
relations: 


(.-a ■ (a. a 


7 

i 

Cp p 

•^n ^ 

(Poje 



+ 



Po' _ 

(7 + 1)%^ 

7-1 

7 + 1 

Po 

J 7 - 1)m/ + 2 


_ 27 m/ - (7 - 1| 


7-1 


(8) 


(9) 


1 % = sin e 


( 10 ) 


The variation of Pq^Pq with the curved portion of the shock will be 

disregarded and the ratio taken as that due to the initial cone. In 
order to detennine the Mach number at any point n on the surface, the 
following equation is used: 



If equation (3) is used. 


Wj^ can be determined and 


(11) 


Hn = sin' 




( 12 ) 


A second condition must be satisfied for the final evaluation of 
the properties of point c . A geometric relation for and is 

obtained by writing the equations for the lines ac and be; thus, 

_ tan(e - b)b - tan ct + yg^ - y^^ 

Xc - 


tan(6 - p )-|3 - tan a 


( 13 ) 
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where 


a 


6 c + 6a 
2 


( 14 ) 


Equation (l 3 ) is solved for several values of a - 6 at the pre- 
ceding point on the body. The results of eqixations (l 3 ) and (5) are 
plotted against where 


<i^a,c = - ^a 


(15) 


The intersection of the two curves represents the point c. The terras p,, 
and Me are obtained by substituting the Intersection values into equa- 
tions (11) and (12). 

In order to obtain a second approximation of the properties of the 
flow at a point on the body, equations (6) and (7) are written; 


<^b,c % - ^b 

^b 


+ w„ 

D C 


0 


c,m 


tan 


l-L 


m- 


^b,m 


®b ®c . ^b 
sin 2 — ^ sin — - 


cos 


^9b + 0c p-b + P(- 


(16) 


(17) 


where the subscript c refers to the intersection values of the first 
approximation. Equation (l 3 ) is rewritten as 


Xb tani 


'■c ,m - 


^6b + Qc ^"b + ^^c^ 


+ ya - 7b - ^a 


tan 


®b ®c 


^'b 




tan Ojji 


Equation (l 4 ) becomes 


(18) 


®c,m “ " ®a 


(19) 
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The results of equations (l8) and (l6) are plotted as against 

dXa^c,m from this plot the intersection values are determined. For 

a body with a negative pressure gradient, the slope at any point on the 
surface must have 0n ^ 9n-l* 

The equations for determining the flow properties of a field point 
and a point on the curved portion of the bow shock by the method of char- 
acteristics have been presented in the literature and are presented for 
the convenience of the reader in appendix A. 

The finite distances between the computed points along the surface 
of the body are determined by the number of q (fig. l) used between 
and q = e on the conical nose. On the afterportion of the body, 

the characteristics are widely separated, the slope of the siurface is 
changing rapidly, and no intersection of the characteristic with the 
surface of the body is obtained behind point t. (See fig. 3(a).) A 
new field point u corresponding to a more closely spaced characteristic 
net (see fig. 3(a)) is obtained by averaging all the properties of the 
points t and r. A more detailed characteristic net used to determine the 
afterportion of the body is presented in figinre 3(b). 

The number of approximations necessary to insure a high degree of 
accuracy was determined by applying the method of this report to a body 
of revolution with zero pressure gradient (a cone) with = 15° and 

Moo = 3-02. The points a, b, d, and so forth, of figure 1 were determined 
from reference 1. The computed 0 (along the surface) rapidly approaches 
its true value asymptotically; that is, a first approximation of the 
point c gave 9^ = lA.82°; a second approximation, 0c, m = 1^.99°; and 
a third approximation, ®c,2m ~ 1^.993°« A second body point was computed 

by using the second-approximation values of points c and e. The results 
were 0j = l4.80° and 0^ m = 14.97°. 

The pressure distribution of a parabolic body used in reference 3 
was determined by the method of cheiracteristics. The resulting pressirre 
distribution was used to determine several ordinates and slopes on the 
initial portion of the nose of the body by the method outlined in this 
report which used the second-approximation values. A comparison of the 
resultant body with the original body is presented in table I and indi- 
cates close agreement. 
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APPLICATION OF THE METHOD TO THE DESIGN OF A SPECIFIC BODY 
Design Specifications and Attainment 


In order to obtain experimental corroboration of the method outlined 
in this report, a model with the following conditions was designed: 

Moo = 5.15, -- = 0.5202, and t] = 15° with an initial cone length 

d(x/L) 

of 1 inch. The x/l, y/L, 0, M, and Cp for the resultant body are 
shown in figure 4. 


The downstream end of the body is defined by the method of char- 
acteristics as the point on the contour where the local static pressure 
becomes zero and can be detemnined from the following relation for a 


constant 


dCi 


d(x/L) ‘ 


L = -a (20) 

— ^ + 1 

n 

d(x/L) 

where Xg, is the ordinate corresponding to the base of the conical 
nose and 


Ap = Px=a - Px=L = Px=a 

The description of the flow by the characteristic net in the 
vicinity of this limiting length is extremely laborious and from the 
nature of the flow has no physical significance. Accordingly, the net 

was terminated at £ = 0.9895 (which will be designated L^) where 
M = 6.205. The maximum radius of the body = O.IOI6 
= 0.6678 and the resultant fineness ratio was 4.919. 



The Bell Telephone Laboratories X-66744 relay computer at the 
Langley Laboratoiy was used for the computations of the example body 
up to 88.8 percent of L. The averaging process previously described 
for determining a more closely spaced characteristic net had to be com- 
puted on desk calculators. This method is very time consuming and 
requires 605 hours for the computation of the entire body (425 hours on 
the Bell machine and I80 hours on the desk calculators). As a comparison. 
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the average time for computing the pressure distribution for a body of 
revolution with a given contour by the method of characteristics on the 
Bell machine is 120 houirs . Convergence on a solution could be expedited 
with a different machine from those available for the computation of 
this project. 

Computation of Cp for Off -Design Mach Nimiber Operation 

In order to determine the effect of varying the free-stream Mach 
number on the body which had a linear pressure distribution at Moo = 5.15, 
the pressure distribution was computed by the method of characteristics 
for Moo = 5.02 and = 4.l6. The ordinates of the body are presented 

in table II. 

The computations for Moo = 5*02 were terminated at the point where 
further computations would have to be done on desk calculators. The 
resultant surface Mach number and pressiore coefficient are shown in 
figures 5 and 6, respectively. The plot of Cp against x/l indicates 

that the pressure distribution is almost linear for Moo = 5.02 but not 
for Moo = 4.l6. A plot of dCp^dx against x/l (fig. 7) is indicative 

of the effect of varying the free-stream Mach numbers. 


In order to illustrate the effect of body shape on the pressure 
distribution, the constant-pressure-gradient model is compared with the 
MCA RM-10 research missile. The contour of the MCA RM-10 model is 
defined as 


where L = 95 inches for comparison with the constant-pressure-gradient 
model. The Mach number and the pressure distribution determined by the 
method of characteristics for the MCA RM-10 research missile at Moo = 5.12 
are shown in figure 8. In figure 9, the slopes of the contours for the 
two bodies are shown where the abscissa is Xy/Lp and is the distance 

from the nose to the point of r^^y . The condition of constant pressure 

gradient requires a relatively blunt afterbody. It is obvious that the 
relation of constant pressure gradient implies an invariant relation 
of 0 with X, and within reasonable limits the value of dp/dx is 
inversely proportional to the length. 


Comparison of Constant-Pressure-Gradient Body With 


MCA RM-10 Research Missile 
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APPARATUS AND TESTS 


Tests on a model of the constant-pressure-gradient body were con- 
ducted in one of the blowdown jets of the Langley Gas Dynamics Branch. 

The test facility consisted of nozzle blocks to give M = 3*05 and 
M = 3.13 with test-section sizes of 8 by 10 inches and 9 "by 9 inches, 
respectively. The test section of each set of nozzle blocks was cali- 
brated by measuring the free-stream static pressiore in the vertical 
center plane with the multitube rake shown in figure 10. High-pressure 
air was heated and discharged throu^ the two-dimensional symmetrical 
nozzle and diffuser to the atmosphere. The pressures were measirred on 
a mercury manometer board and photographically recorded. The resilLts 
shown in figures 11 and 12 indicate that M = 3-05 t0.045 and M = 5«13 
±0.0^9, respectively, over the entire region that was surveyed. A 
schlieren photograph of the flow in a region of the test section of the 

M = 5.05 nozzle at a Reynolds number per inch of 2.77 X 10^ is shown 
in figure 13 and a schlieren of the flow in the M = 3-13 nozzle at a 

Reynolds number per inch of 2.6k X 10^, in figure l4. 

Measurement of aerodynamic characteristics that are affected by 
the boundary -layer development are governed by the transitional Reynolds 
number. The transitional Reynolds nimiber was obtained by heat-transfer 
measurements on a 10^ cone to indicate the steadiness of the free stream 
in this facility as compared with that of other NACA facilities. The 
method and results are discussed in appendix B. 

The constant-pressure-gradient model has a length of 17* 55 inches 
= 0.9781^ and a maximum diameter of 3.608 inches and is constructed 

of Paraplex with a steel central core and coated with heat-resistant 
phenolic lacquer for added strength and heat resistance. A sketch of 
the model and a table of ordinates are shown in figure I5 and a photograph 
of the model in figure 16. The completed model was accurately measured 
on a jig-boring machine and the resultant ordinates are plotted as a 
deviation of slope of the surface from the slope of the prescribed ordi- 
nates in figure I7. The differences between the prescribed and measured 
slopes are shown to be small. The model was mounted on a sting which 
was supported by a strut spanning the tunnel downstream of the diffuser. 

A sketch of the model and sting support momted in the M = 3.13 timnel 
is shown in figure I8. Thirty-five static -pressure orifices ( 0.020-inch 

inside diameter) were located at ^inch increments along a longitudinal 

element of the body with one orifice diametrically opposite the foremost 
orifice to determine the angle of attack. The pressures were measured 
on a multitube mercmy manometer board and photographically recorded. 
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The effect of orifice size on the static -pres sure measurements on 
a body of revolution vith a small boundary layer was determined on a small 
50*^ cone. Three orifices (0.004-^ 0.0135-> ond 0.0200-inch inside diam- 
eter^ including the orifice size for the design model) were located at 
two axial stations corresponding to the locations of the first two static - 
pressure orifices on the constant-pressure-gradient model. A sketch of 
the model illustrating the size and location of the orifices is shown 
in figure 19 and a photograph of the model, in figure 20. The cone was 
tested with the orifices at 0 = 0^ and 90^ at M^o = 5.15 over the 
same range of free-stream Reynolds numbers used for the tests of the 
design model. A comparison of the measiored Mach number with the predicted 
Mach number (ref. l) indicated that the error is a maximum for the small- 
est orifice and is essentially zero for the orifice with a 0.020-inch 
inside diameter. 

The effect of free-stream static -pres sure variations on the 
longitudinal -pressure distribution for the design model was determined 
by rotating the model about its axis of symmetry and recording the pres- 
sures for 0=0^, 90^, and l80^. In the M = 3.05 tunnel, the model 
was also tested at the three axial stations shown in figure 11 with the 
orifices at 0 = 0^ to determine fiorther the effect of flow discontinu- 
ities. The ranges of test Reynolds numbers based on body length were 

3.61 X 10'^ to 4.86 X 10'^ for M = 3.05 and 3.76 X lo"^ to 4.62 X 10^ 
for M = 3.13. Composite schlieren photographs taken dioring the tests 
are shown in figures 21 and 22. 


ACCURACY OF TEST DATA 


A number of pressures were measured for the purpose of determining 
the free-stream Mach number and the pressiure coefficient along the model. 
The overall accuracy is dependent upon the individual accuracy of each 
measurement. The accuracy of measuring the stagnation pressvire on a gage 
was ±0.1 lb/in.2, and the accuracy of reading it on a manometer board was 
±0.1 inch of mercury. The overall accuracy for determining Cp is ±0.003. 

A comparison of the experimental points with the computed curve indicates 
a deviation greater than the accuracy of measuring Cp. This deviation is 

explained by the deviation of actual model contour from the design (see 
fig. 17) and the growth of boundary layer. 


RESULTS AND DISCUSSION 
Tests at M = 3.05 


The schlieren photograph (fig. I 3 ) shows several compression waves 
symmetrical about the horizontal center plane of the tunnel. The nozzle 
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calibration indicates that only one set of these waves has sufficient 
strength to be of concern. The model was tested in three axial positions 
(fig. 11) and three radial positions with the model in position 1. The 
results of these tests are plotted as Cp against x/l in figure 25 

where the solid line is the computed distribution for Moo = 3.02. The 

region of compression causes a pressure rise at = 0.195 to 0.31 with 

the model in position 1 and is moved forward to ^ = O.O 85 to 0.20 with 

the model in position 2 in the tunnel. (These bumps are dotted for 

clarity of fig. 25(a).) In axial position 3 the compression does not 
intersect the body. The result of rotating the body about its axis of 
symmetry in position 1 is shown in figirre 23(b). The pressure rise that 
occurs when the line of orifices intersect the compression wave is dotted 
for clarity. The slope of the compression can be determined from the 
schlieren (fig. 13) and cam be checked with the results of rotating the 

model about its aixis. The Reynolds number was varied from 56.1 X 10^ 

to ^4-8. 6 X 10^ with the model in position 1 and the orifice at 0 = 0°. 
There was no effect on the pressure distribution (fig. 25(c)). 

The experimental data (fig. 23) were corrected for the free-stream 
static -pressure variation by simple superposition except in the region 
of the bump. The bump in the Cp curve due to the compression waves in 

the nozzle was corrected by decreasing the measirred static pressirre by 
the ratio of the pressmre rise across an oblique shock, the intensity of 
which corresponds to the decrease in Mach nimiber shown in figure 11. 

The coarseness of the nozzle calibration allows considerable latitude in 
the detail of the Mach number distribution in the compression region and 
the distribution was assumed so that it was compatible with the calibra- 
tion and eliminated the bimip in the Cp curve. The effect of this com- 
pression had been completely eliminated when the compression distribution 
in the flow was applied as a correction to the model in all positions, as 
is shown in figure 24. An inspection of figure 24 shows an excellent 
agreement between the computed and experimental Cp curves. 

The pressure rise in the vicinity of the base of the body is asso- 
ciated with the boundary-layer separation as seen in figure 21. 


Tests at M = 5.13 

The body was rotated about its axis of symmetry to evaluate the 
effect of the free-stream discontinuities on the longitudinal -pressure 
distribution. The experimental data corrected for the free-stream static - 
pressure variations by simple superposition are shown in figure 25(a) with 
the computed pressinre distribution. The first two orifices located dia- 
metrically opposite from each other, 0.5 inch from the nose, have a 
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reversal of pressure differential for (f) - qP and l80^; this reversal 
indicates that the model was at a slight angle of attack with the line 
of orifices on the windward side for 0 = 0^. The Reynolds number was 

varied from 37*6 X 10^ to 46.2 X 10^ with the orifices at ^ = 0^ with 
no effect on the pressure distribution (fig. 25(b)). 

A close examination of the flow schlieren (fig. 22) shows the 
conical shock angle is greater than the computed value from reference 1^ 
but the accuracy of measiorement is not sufficient to determine the free- 
stream Mach mmiber. A comparison of the boundary-layer thickness along 
the nozzle surface with that in figure l4 shows the boundary layer is 
approximately 0.032 inch thicker when the model is mounted in the tunnel. 
The effect of the model in the small test section might be to reduce the 
free-stream Mach number from 5.13 to 3.07. When M = 5.07 is used to 
reduce the experimental data^ there is good agreement with the computed 
pressure distribution for the design conditions. The pressure rise in 
the vicinity of the base of the body is due to the boimdary-layer 
separation. 


CONCLUDING REMARKS 


A method based upon the theory of characteristics has been devel- 
oped to compute the contour of a body for a prescribed pressure distri- 
bution. An indication of the acciiracy of the method was obtained by 
using the pressure distribution which had been determined by the conven- 
tional method of characteristics to compute the ordinates near the nose 
of a body. The agreement between the original and calculated ordinates 
was good. The contoirr of a body with a constant pressure gradient at 
a Mach number of 3.13 was computed; a model was constructed and tested 
at Mach mmibers of 5.05 and 3.13. Excellent agreement with the computed 
distribution was obtained. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 21, 1955- 
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APPENDIX A 

EQUATION FOR DETERMINING THE FLOW PROPERTIES OF A FIELD POINT 

The equations for determining the flow properties of a field point 
nnf^ a point on the curved portion of the bow shock have been given in 
many references and are presented in this appendix for the convenience 
of the reader. 

The equations for the first approximation of a field point e (fig. l) 
ane as follows: 


yc - yd - ^c ^^(^c ^^c) + ^d ^^(^d - ^d) 

tan(e^ - Pd) - tan(0(, + pj 


(Al) 


7e = (Xe - XeL)tan(0(i - p^) + 


(A2) 


rri|r3 V r* / 

1 - ^ - (»c - I'd + ^ (’■e - =‘d) - ^ 

— tan p„ + tan p , 

% ^ ^ 


) ^ 
-) W. 


where 


sin [X tan sin 0 
cos (9 + [i) 


(A5) 


(A4) 


Then, 


Wg = W(, + Wg 


(0e - ®c)^^ ^ (^e 



(A5) 


The term Pg is determined from equation (l2). 

The equations for the second approximation are as follows: 
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yc 




+ e. 




^e,m 


tan 


’6d ■*■ ®e l^d ■*■ 


2 / 


- tan 


/6c + 0e Me + ^^e^ 

[—r-^—T—J 


(A6) 


ye,m = (xe,m - xa)tan ^ ^ ^ 


(A7) 


. (9e - 0d)tan 

W<i + Wg ° 2 


, /.. ,.,x ™d,m 2(Wg + Wg) x^ - Xg 

^ ■ IT^' ; w; - V, 


yc + ye 


d0^ 


tan + tan 

2 Wd ^ ^e 2 


(A8) 


sin 


^^r. + M, 


tan 


sin 


0+0 
_£ e 


c ,m 


/®c + ®e 

COS + 


(A9) 


®e,m *^^ 0 , 6 , in ■*■ ®c 


(AlO) 


W, 


e.m 


(9e,m - 0 c) tan 




yc + Ye 
2 


W + W 

^ o ^ + W„ (All) 


The final x-value of the field point is that calculated from equa- 
tion (a 6) where Xg^j^ corresponds to an average value of p and 0 for 

the first approximation obtained by using points c, d, and e. 

In order to determine the properties of point h, the intersection 
of a Aj characteristic with the bow shock, assume the line fh (fig. l) 

is a straight line, the slope of which is tan e, and the line gh is 
straight with a slope of tan(p + 6)g. As a result. 
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Xh = 


'8 


- + Xf tan e - x„ tein(n + 6) 


'g 


tan € - tan(|i + 0) 


(A12) 


and 


= (""h - Xf)tan € + Yf 


(A15) 


With the use of the general equation of a shock. 


W 


6>g 


W, 


1 + 


g 


Jeg Wg 


dGg^h tan Ug - Ig = 0 


(Aik) 


Then 


d0 


8,h 


Xh - Xg 

yg 


W. 


Zg+ 1 


^^8 


W, 


g 


ds/e Wg 


tan 




(A15) 


where the subscript 9g corresponds to a deviation across the shock 
at h equal to 0 of point g. Therefore, 


®h = ^h = ®g + <^®g,h 


(Al6) 


and 


Wv, = Wa „ + (M) d0 


^ ll5 


’g,h 


(A17) 


g 


The method of determining (dW/d6)g and will be presented 

after the equations for the second approximation of point h which are 
as follows: 


"h. 


Yg - yf + Xf 


tan 


e + €v 


- x_ tan 
& 




m 


^ ^h ^ ^^g + ®g + M-h + ®h 

- tan — 


(A18) 


tan 


2 


2 
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(A19) 



g 


(A20) 



(A21) 


The term is computed from equation (A9) hy substituting g for c 

and h for e. The terms (4^1 and Wq _ are constant for both 

\d6/0g 

first and second approximations. The shock angle corresponding to 

is used to determine the next point on the shock but not to determine 
the location of point h^. 

As outlined in reference 2, dW/d6 is the tangent, at the point 
6 = 9g, to the curve that represents W as a function of 6. The usual 

procedure is to determine this property graphically and regard the quan- 
tity as constant along the curved portion of the shock. An analytical 
expression is necessary to determine dW/dS to a more accurate degree 
and to consider the curvature of the bow shock for computing the char- 
acteristic net on a Bell Laboratories X-66744 relay computer. If the 
two-dimensional shock equation is used, 



(A22) 


and 



Since 



then, by using the notation of this report 
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d 1 _ / y + 1 ^^£2 A 1 

de tan 5g 2 _ i j cos\ 


2 . 


4 - 1 
2 




2 sin^e 


g 


1 

de _ dS 

<aSg " A. 1 

de tan 6g 


(A25) 


(A26) 



sin 




(y + 1 ^an e^ 

2 tan (Eg - 6g) 



(A27) 
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The variation of dW/d5 from the point where shock curvature begins to 
the last point on the shock where a characteristic of Ajj family inter- 
sects the stern of the body (see fig. 26) was computed for the constant- 
pressure -gradient model. For Moo = 5-15> dw/d6 varied from -0.55957 

at V, to -0.511^5 at V T . 

1 n+1 
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APPENDIX B 


METHOD OF DETERMINING THE TRANSITIONAL REYNOLDS NUMBER 
OF THE M = 5.13 NOZZLE 


In a project of this type^ where a computed distribution is to be 
compared with the actual pressure distribution on a body of revolution^ 
it is interesting to know the nature of the boundary layer because of 
its effect on the model contour. 

The boundary-layer transition Reynolds number on a 10^ cone has 
been used as a measure of the steadiness of the free stream in various 
NACA facilities (ref. 4). A similar highly polished thin-walled cone 
was used to determine the R^^ by surface-temperature measurements in 

the M = 5*15 jet. The high free-stream Reynolds number per foot pro- 
duces an extremely thin boimdary layer which requires that the ratio of 
average height of surface roughness to height of boundary layer be small. 
The ratio is made small by a high degree of polish on the model surface 
and by reducing the free-stream Reynolds number. 

The model with a 0.050-inch skin thickness was machined from a 
single piece of stainless steel. The surface was successively ground 
with No. 1 / 0 , No. 2 / 0 , and No. 3/0 fine emery papers and the direction 
of grinding with each paper was alternated from circumferential to axial^ 
and conversely. A cutting oil of 5 percent paraffin with filtered kerosene 
was used with each grinding. The final polish was obtained with a com- 
mercial polish which contains a chromium grit whose average size is 
1 micron. 

The temperature distribution along the surface of the model was 
determined by I 8 iron-constantan thermocouples installed in a 20 ^ helix 
about the conical body with the first thermocouple located 5 A i^ch from 
the nose as shown in figure 27 and a picture of the completed model in 
figure 28 . A continuous record of the s\irf ace-temperature variation was 
made throughout the test with a multi-channel recording oscillograph. 

In order to facilitate a closer inspection of the boundary layer 
by visual means ^ a schlieren system similar to that in reference 5^ which 
permits magnification greater in a plane normal to the flow than in the 
direction of the flow^ was constructed. The system incorporates three 
spectacle-size lenses: one spherical and two cylindrical. The optical 

diagrams are shown in figure 29- A magnification of 6^ normal to the 
flow^ was obtained with this lens system incorporating a +^k diopter 
cylindrical lens with a horizontal axis (lens no. 4 in fig. 29 ) and a 
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+6 diopter cylinder with a vertical axis (lens no. 5 in iig- 29). The 
lens mount in figure 50 allows the lens position and cylinder axes to be 
varied. A comparison of the 10° cone using a conventional schlieren with 
this system is shown in figure 31. A 1-inch grid is superimposed on the 
asymmetrical schlieren. Since the lens system was not achromatic^ a 
filter which transmits the yellow portion of the spectrum was used. 

A second minimum was installed to reduce the stagnation pressure and 
thereby reduce the free-stream Reynolds number Rco and move the tran- 
sition point onto the instrxmiented portion of the cone and to a point 
where the boundary layer has a substantial thickness. The diffuser con- 
sisted of two 11° wedges with a length adjusted to avoid choking. The 
portion of the diffuser downstream of the second minimimi was extended to 
a point beyond the model support as shown in figure 32. This extension 
allowed the stagnation pressure to be reduced from l60 Ib/in.^ to 
65 lb/in.2 gage and a corresponding reduction of Reynolds number per inch 
from 2.15 X 10^ to 0.894 X 10^. 

The variation of the transitional Reynolds number with free-stream 
Reynolds mmiber is shown in figure 55- The transition point corresponding 

to this curve is accurate to inch. 

The transition point'was also determined from the flow schlieren 
(fig. 5^) aJid the corresponding transitional Reynolds number is plotted 
in figure 55- The tests indicate that the transition point on the 
constant-pressure-gradient model is so far forward that the boundary 
layer is turbulent on the entire model. 
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TABLE I.- COMPARISON OF ORDINATES CCMPUTED FROM THE KNOWN PRESSURE 
DISTRIBUTION WITH THE ACTUAL ORDINATES OF A 
PARABOLIC BODY AT Mco = 3-02 


Method 

in. 

y , in. 

0, deg 

Pressure 

1.0354 

0.55618 

17.540 

Computedl 

1.0354 

.33620 

17.541 

Pressure 

1.0785 

.34979 

17.500 

Computedl 

1.0785 

.54981 

17.503 

Pressure 

1.1256 

.36400 

17.462 

Computed^ 

1.1236 

.56401 

17.464 


^Con 5 >uted by the method of characteristics with 
contour of the body defined by the equation 



where 


L = kO. 
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TABLE II.- ORDINATES FOR THE CONSTANT-PRESSURE-GRADIENT MODEL 
DETERMINED BY THE METHOD OF CHARACTERISTICS 


x/l 

r/L 

0.015900 

0.0057244 

.014495 

.0058858 

.015120 

.0040507 

.015777 

.0042262 

.016472 

.0044115 

.017211 

.0046081 

.017999 

.0048175 

.018845 

.0050419 

.019754 

.0052853 

.020740 

.0055447 

. 021815 

.0058292 

. 022995 

.0061412 

.024301 

.0064860 

.025760 

.0068705 

. 027405 

.0075054 

.029286 

.0077975 

.051470 

.0085700 

. 034058 

.0090465 

.037204 

. 0098667 

.041172 

. 010897 

.042472 

.011234 

. 044254 

.011689 

.046251 

.012210 

. 048574 

.012757 

.050616 

.013335 

.053003 

.013944 

.055549 

.014595 

.058285 

.015291 

.061241 

.016042 

.064445 

.016855 

.067957 

.017738 

.071825 

.018708 

.076121 

.019781 

.080945 

.020980 

. 086427 

.022532 

. 092729 

.023877 


x/l 

r/L 

0.10012 

0 . 025672 

. 10896 

.027797 

. 11986 

.050383 

.15585 

.053645 

. 14475 

.036155 

.15204 

.037806 

.15979 

.059545 

. 16805 

.041379 

. 17696 

.043330 

.18659 

.045408 

. 19707 

.047637 

. 20858 

.050041 

. 22126 

.052638 

. 25542 

.055471 

.25155 

.058576 

.26949 

. 062000 

.29042 

.065802 

.51500 

.070058 

. 54442 

.074843 

.38064 

. 080258 

.42679 

.086351 

.48878 

.093003 

•57912 

.099151 

.61182 

. 10025 

. 66079 

. 10058 

.72415 

.098326 

.80425 

.090054 

.85734 

. 080054 

.88806 

.072029 

.92456 

.059352 

.95000 

.047426 

.96611 

.037746 

.97610 

.030486 

.98209 

.025457 

. 98742 

.020324 

.98954 

.018061 


h 



MCA TN 3555 


ro 

ON 



Figure 2.- Assumed axial variation of pressure coefficient for body of revolution. 


MCA TN 3555 





(a) Station 44 to 71 . 1 . 

Figure 3*- Characteristic net for afterportion of body. 
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Figure 5*- Concluded. 
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Figure 4,- Axial variation or Cp^ and 0 for constant -pressure-gradient body at 

Meo = 3.13. 
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Figure 5*- Axial variation of surface Mach number for a body designed to 
have constant dp/dx at = 3.13 when operating at M^^, = 3.02 
and Mqo = 4.l6. 
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Figure 6.- Axial variation of pressure coefficient for body of figure 4 
when operated at = 5.02 and = 4.l6. 
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Figure Variation of dCp/dx for linear constant -pressure-gradient body at designed condition 
of Moo = and off-design condition of Moo = 3-02 and Moo = ^.l6. 
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Figure 8.- Distribution of M and Cp for MCA RM-10 research missile 

at = 5.12. 
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Figure 9 .- Comparison of slope of surface of the NACA RM-10 missile and the constant-pressure- 

gradient model. 
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Figure 11.- Nozzle calibration for M = 5-05 at a Reynolds number per inch of 2.77 x 10 . 
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Figure 12.- Nozzle calibration for M = 5*15 at a Reynolds number per inch of 2.6k x 10^. 
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Figure IJ.- Schlieren photograph of M = 5.O5 test section. 
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Figure Ik.- Composite schlieren photograph of M = 3.I5 test section. 
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Figure 15.- Sketch of constant -pressure-gradient model with pertinent dimensions. 
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Figure I 7 .- Design slope for the contour of the model and deviation of the experimental model 

from design. 
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Figure l8.- Sketch of model in M = 5*15 tunnel. 
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Figure 19*- Sketch of cone used in orifice investigation. 
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Figure 20.- Photograph of cone used in orifice investigation. 
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Figure 21.- Composite schlieren photographs of constant-pressure-gradient model at M = 5-05 

for Reynolds nimibers of 0 and 48.6 x 10^. 
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Figure 22.- Composite schlieren photograph of model at M = 3-15 for Reynolds numbers of 0 

and 46.2 X 10^. 
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(a) Effect of varying axial position of model in tvinnel 
at Reynolds number 48.6 X 10^ with orifice at 0°. 

Figure 25.- Axial distribution of pressure coefficient of constant -pres sure - 
gradient body at Moo = 5-05. Flagged symbols denote orifice displaced 
l80°. Dashed lines indicate the effect of the compression. 
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(b) 


Effect of radial position of the orifices with the model at position 1 
in tunnel at Reynolds number 48.6 x 10^. 


Figure 25.- Continued. 
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(c) 


Effect of varying Reynolds number with the orifices at 


0° at position 1. 


Figvire 23-- Concluded 
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(a) Effect of varying axial position of model in t\mnel 

at Reynolds number 48.6 x 10^ corrected for free-stream 
disturbances . 

Figure 24.- Axial distribution of pressure coefficient of constant-pressiire- 
gradient body at = 5*05 corrected for free-stream disturbances. 

Flagged symbols denote orifice displaced 180°. 
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(b) Effect of radial position of the orifices with the model at position 1 
in tunnel at Reynolds number 48.6 x 10^ corrected for free-stream 
disturbances . 


Figure 24.- Continued. 
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(c) Effect of varying Reynolds number with the orifices at 0° at posi- 
tion 1 corrected for free-stream disturbances. 


Figure 2k.- Concluded. 
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(a) Effect of radial position of orifices at Reynolds 
number 46.2 X 10^. 

Figure 25.- Axial distribution of pressure coefficient of constant -pres sure - 
gradient body at = 3*15 corrected for free-stream disturbances. 
Flagged symbols denote orifice displaced l80^. 
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(b) Effect of varying Reynolds number with orifices at 0°. 
Figure 25 .- Concluded. 
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Figiare 27*- Sketch of the model used for the transition investigation. 
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Figure 28.- Photograph of 10° cone used for transition investigation. 
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Vertical view 


Figure 29«- Optical lens system for producing greater magnification in a plane normal to the flow. 
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(a) Asymmetrical. 



(b) Conventional. 
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Figure 51 •- Comparison of 10° cone using conventional and asymmetrical 

schlierens. 
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Figure 52.- Sketch of model and diffuser for transition tests. 
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Figure 53.- Variation of transitional Reynolds number with stream Reynolds number measured with 

a 10^ cone . 
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(a) Asymmetrical. 



(b) Conventional. L-9051? 

Figure 54.- Con^arison of conventional and asymmetrical schlierens on a 

10° cone at R = 12 x 10^. 
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